arXiV:math-ph/0503074vl 31Mar 2005 


ON THE COMPLEXITY OF SOME BIRATIONAL 
TRANSFORMATIONS. 

J. CH. ANGLES D’AURIACf, J. M. MAILLARD| AND C. M. VIALLET* 


Abstract. Using three different approaches, we analyze the complexity of 
various birational maps constructed from simple operations (inversions) on 
square matrices of arbitrary size. The first approach consists in the study 
of the images of lines, and relies mainly on univariate polynomial algebra, 
the second approach is a singularity analysis, and the third method is more 
numerical, using integer arithmetics. Each method has its own domain of 
application, but they give corroborating results, and lead us to a conjecture 
on the complexity of a class of maps constructed from matrix inversions. 


1. Presentation 

The investigation of birational representations of Coxeter groups acting on pro¬ 
jective spaces of various dimensions appeared some years ago to be of interest to 
understand the structure of lattice models of statistical mechanics mm Birational 
dynamical systems have also been studied for their own sake with various methods 
ranging from analysis to algebra. A common ingredient to these subjects is the 
study of iterations of infinite order birational transformations, and in particular 
their complexity, measured by the rate of growth of the degree of their iterates; see 

for example IS1I1I11E1IZ1IH1I3IHIII1II1II1 

We perform this analysis for a definite class of transformations, defined from 
elementary operations on matrices of size qxq, the entries of the matrices being the 
natural coordinates of complex projective spaces CP„. Depending on the specific 
form of the matrices, the dimension n will take different values {n < q^ — 1). 

We explain, exemplify, and confront three different approaches to the problem. 
We also present a conjecture for the value of the complexity for a family of trans¬ 
formations of interest to statistical mechanics. 

The paper is organized as follows. We state in section[2lthe problem of calculating 
the complexity of a birational transformation acting on a projective space, and 
define the basic objects of interest. We introduce four families of maps, which 
will be used for explicit calculations. In section El we indicate how to infer the 
generating function of the sequence of degree of iterates of a map from its first 
terms. This provides a first method of calculation of the complexity. In section^ we 
calculate exactly the sequence of degrees by an analysis of the singularity structure 
for one of the families of maps. In sectional we describe an arithmetic approach, 
where we examine the action of iterates on rational points (integer homogeneous 
coordinates), and simply measure the growth of the size of the coordinates. This 
yields approximate values of the complexity. We conclude with a conjecture. 
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2. The problem 

Let isT be a birational transformation of complex projective space CP„. If we 
write K in terms of homogeneous coordinates, it appears as a polynomial transfor¬ 
mation given by n -I- 1 homogeneous polynomials of the same degree d. With the 
rule that we should factorize out any common factor, d is well defined in a given 
system of coordinates. Of course it is not invariant by changes of coordinates. We 
may construct the sequence {dn}, of the degrees of the iterates Lf" of K. 

The growth of the sequence is a measure of the complexity of K. In the 
absence of factorizations of the polynomials the sequence would just be 

(2.1) dn = d^ = d”. 


What happens is that if some factorizations appear, they induce a drop of the 
degree, so that we only have an upper bound 

(2.2) d„ < d”. 

The drop may even be so important that the growth of d„ becomes polynomial 
and not exponential anymore. A measure of the growth is the algebraic entropy 

(2.3) e= lim ilogd„, 

n —Ti 

or the complexity 

(2.4) A = exp(e). 


Both the entropy e and the complexity A are invariant by any birational change 
of coordinates. They are canonically associated to the map K. Our aim is to 
calculate them for definite classes of maps, which we now describe. 

Suppose M is a. q X q matrix, and consider the two simple rational involutions 
/ and J: the involution / is the matrix inverse up to a factor (i.e. when written 
polynomially it amounts to replacing each entry by its cofactor). The involution 
J is the element by element inverse (also called Hadamard inverse, which replaces 
each entry My by its inverse 1/My). The two involutions / and J do not commute, 
and their composition K = I o J is generically of infinite order. 

The map K acts naturally on <CPq 2 _i. It is however possible to define various 
reductions to smaller projective spaces in the following wa,v[Tl|. For a given size 
of square matrices, we define a pattern as a set of equalities between entries of 
the matrix. The set of all pattern is the set of all partitions of the entries of the 
matrix. An example of a pattern is “all diagonal entries equal, all off-diagonal entries 
equal”. This corresponds to the partition of the entries in two parts (diagonal + off- 
diagonal). Clearly any pattern is preserved by the action of J. We call admissible 
a pattern which is also stable by / (or equivalently K). 

All admissible patterns have been classified for 9 = 4 and some of them for 
g = 5 in It has been also shown that A can vary considerably from 

one admissible pattern to another. For example for 5x5 cyclic and symmetric 
matrices one has A = 1 (polynomial growth), whereas with the cyclic matrices one 
gets A = (7-b 375)72. 

We will focus on four fundamental admissible patterns, which exist whatever the 
size g X g of the matrix is. The first one is the pattern {S) of symmetric matrices. 
The second one (C) is the pattern of the cyclic matrices defined by Mij = M^+i 
(with indices taken modulo g). The third one is the pattern of matrices which are 
at the same time cyclic and symmetric (CS). The last one is the general pattern 
(G), without equality conditions between the entries. 
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3. A FIRST APPROACH: GENERATING FUNCTIONS 

From the sequence of degrees {dn}, it is possible to construct a generating func¬ 
tion 

CXD 

(3.1) /(u) = ^d„M". 

n—0 

Since the degrees are bounded by (E3, the series isu always has a non zero 
radius of convergence p. Actually 

(3.2) p=j, 

The calculation method is the following: calculate explicitly the first terms of 
the series, and try to infer the values of the generating function. The method is 
sensible if the generating function is rational. 

The striking fact is that indeed the generating function f{u) happens to be a 
rational fraction with integer coefficients in most cases. The consequence is that a 
finite number of terms of the series determine it completely. For reversible maps 
(i.e. when there exists a similarity relation between the map and its inverse), we 
have not found any counterexample to this rule. There are however non-reversible 
maps for which the generating function is not rational m Another consequence of 
the rationality of / is that A is an algebraic integer, and we have no counterexample 
yet to that. 

For practical purposes, it is necessary to push the calculation of the degree of 
the iterates as far as possible. Instead of evaluating the full iterate, it is sufficient 
to consider the image of a generic line I with running point 

(3.3) l{t) = [ao + bo t, ai + bi t,..., Qn + bn t], 

where ai,bi are arbitrary coefficients, and evaluate the images of l(t) by if”. The 
degree is read off from this image. The calculation may furthermore be improved 
by using integer coefficients in and calculating (formal calculation software are 
quite efficient at that) over polynomial with coefficients in Z/Zp with p a sufficiently 
large prime integer. Taking different values of p and of the coefficients ai,bi helps 
eliminating the accidental simplifications which may occur. 

Suppose we have the degree dn for the first values of n, say n = 1... Ujuax- We 
may fit the series with a Fade approximant F, with numerator (resp. denominator) 
of degree N ( resp. M), such that 

(3.4) N + M = rimax ~ 1 

N running from 0 to n„,ax ~ 1- Our experience is that, if Umax is large enough, 
the rational fraction F we find simplifies drastically, and stabilizes for some cen¬ 
tral values of N. This usually means that the exact generating function has been 
reached. 

Note that the expansion of the non optimal [N, M] Fade approximants yield non 
integer, or negative coefficients in the expansion of F, in contradiction with these 
coefficients being a degree. Table Q] displays the “exact” expression we have inferred 
for the generating function for various values of q for the (CS) pattern, as well as 
the value of to = -I- M and the value of rimax- 

When Umax is larger than to, we have a prediction on the next values of the 
degree, and this gives confidence that the result is exact. 

In Table QJ we also give the inverse of the modulus of the smallest zero of the 
denominator, as well as a numerical value computed as explained in section 0 
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q 

fq{u) 

^max 

m 

A 

'^num 

4W 

(l+uf 

{1-uf 

OO 

4 

1 


5W 

(^l+u+2u'^'j 

14 

9 

1 

1.0062 

{l—uf {1+u+u'^) 

6 

{l+2uf 

15 

4 

4 

4.0003 

(1—u)(l—4u) 

7W 

(^l-|-u-|-3u^) 

12 

9 

6.854102 

6.8541 

{l—u){l+u+u'^){l—7u+u'^) 

8 

{l+u)(^l+2u—u'^^ 

11 

7 

10.331852 

10.3317 

{l—u){l—llu+7u^—u^) 

9 

(^l+u+3u‘^—3u^^ 

11 

13 

12.832689 

12.8326 

{l—u){l—l?>u+2vP‘+u^+\2u'^—'S>u^+u^) 

10 


9 

5 

17.944273 

17.9453 

(1—u)(l—IS m-I-u^) 

iiW 

(^1-|-u-|-5m^) 

7 

9 

22.956439 

22.9562 

( 1 —It) (1) (1—23m-|-u^) 

12 

(^l-|-4u— 3it^) {\+2u—v?^ 

8 

8 

25.812541 

25.8105 

(l-u)(l-27u-f3lM2-9u3) 

13(*) 

(^l-|-u-|-6u^) 


9 

33.970562 

33.9719 

( 1 —rt) (1) (1—34m-|-u^ ) 


Table 1. Generating functions for the cyclic symmetric {CS) pat¬ 
terns. The formulae for prime values of q, tagged with a (*), can 
be proved. Umax is the maximum number of iteration performed, 
m refers to the Fade approximation, A is the complexity, Anum is 
the numerical complexity calculated in section 0 
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4. A SECOND APPROACH: SINGULARITY ANALYSIS 

In this section we prove that the complexity of the patterns (CS) for prime q is 
a quadratic integer, by showing that the sequence of degrees verifies a linear recur¬ 
rence relation of length 2 with integer coefficients. This implies that the generating 
function of the degrees is a rational fraction and corroborates a part of the results 
given in tabled 

4.1. Some notations. Let M be a cyclic symmetric matrix of size q x q. The 
matrix M may be written in terms of the basic cycle of order q: 


/ 

0 

1 

0 ••• 

0 



0 

0 

1 • • • 

0 



O ; 

O ; 

O ; 

1 


V 

1 

0 

0 ••• 

0 

/ 


A^even = Xq + Xi {a + H-h Xp-i p=^ + l 

Modd = Xo + Xi {a + H-h Xp-i ^ ^ ^ 9^ 

when q is even and odd respectively. 

The parameter space is a projective space CPp_i of dimension p — 1, with p = 
q/2 -b 1 if g is even and p = {q + l)/2 ii q is odd. We use homogeneous coordinates 

[xq, ■ ■ • , Xp_i]. 

We will study the two elementary transformations / and J acting on M. Both 
are rational involutions (and are thus birational transformations). 

The Hadamard inverse J may be written polynomially in terms of the homoge¬ 
neous coordinates 

(4.1) J : [xo,...,Xp_i]—> Xk] 

k^l k^p—1 

The matrix inverse /, up to a factor, transforms cyclic matrices into cyclic matrices, 
and symmetric matrices into symmetric matrices. It thus acts on cyclic symmetric 
matrices. 

For cyclic symmetric matrices, the matrix inverse I is related to the Hadamard 
inverse J, by a similarity transformation: 

(4.2) I = C-^oJoC 

The transformation C acts linearly on the p homogeneous coordinates. Denoting w 
the g-th root of unity, C is given by the px p matrix with entries: 

(4.3) a.o = i, a.. = + ^^0 

for q odd and 

(4.4) a.o = 1, a.. = a;’'* + 4l r^0,p-l, 

Cr,p-1 = (- 1 )" 

for q even. 

The matrix C verifies = 1 
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4.2. Sequences of surfaces, sequences of degrees. Consider now a sequence 
of hypersurfaces in CPp_i, obtained by applying successively I, then J, then I, and 
so on, starting with a generic hyperplane Sq. Each surface S'„ has a polynomial 
equation, of degree (i„, which we also denote S'„. Since for non singular points, 

(4.5) X € S2n ^ J{x)€S2n-l, 

S 2 n can be obtained from S 2 n-i by substituting the coordinates of x with the 
homogeneous polynomial expression of the coordinates of J{x) in S' 2 Tt-i(a:)- Notice 
that, since J is an involution, S 2 n-i may be obtained from S 2 n in the same manner. 

What happens at the level of the equations is that S 2 n-iiJix)) may factorize. 
One of the factors is S 2 n{x). The only other possible factors are powers of the 
coordinates of a; = (xq, • • • , Xp-i) as explained in the lemma below. Relation 

P-i (fc) 

(4.6) S' 2 „-i(J(a;)) = S 2 „(x) ■ 

k=0 

defines the exponents a^n-i. 

4.3. A lemma. The previous relation is crucial. Its proof is elementary and goes 
as follows. 

Suppose B is a birational involution. When written in terms of the homogeneous 
coordinates, B'^ appears as the multiplication by some common polynomial factor 
of all the coordinates, that is to say the identity transformation in projective space. 

(4.7) B{B{x)) = [Ksix) ■ xq, kb{x) • 0 : 2 , • ■ ■, kb{x) ■ Xp-i] 
with kb{x) = some polynomial. 

We then have, if two algebraic hypersurfaces S and S’ are the proper images of 
each other by involution B-. 

(4.8) S{B{x)) = S'{x)-R{x) 

(4.9) S'(B{x)) = S{x)-T{x) 

with R and T some polynomial expression of the coordinates. We then have us¬ 
ing 114.811 and 114.911 : 

(4.10) KB(a:)'^®s(S) . s{x) = S{x) ■ R{B{x)) ■ T{x) 
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that is to say 

(4.11) = R{B{x)) ■ T{x) 

Equation 114.1 111 shows that the only factors in the right hand side of equation II4.9II 
are the equation of S', and polynomial expressions T{x) which divides kb{x), pos¬ 
sibly raised to some power. 

In the specific example B = J, and S = S 2 n-i, using 



p-1 

(4.12) 

i^jix) = n 



we get : 

p-1 

(4.13) 

S2n-l{J(.x))=S2n{x)-l[ 
2 = 0 


with Xi{x) the i’th coordinate of t, and pi some integer power. 
This ends the proof of formula (ESI). 


4.4. Recurrence relation. Similarly to equation E31, we have: 

p-i 

(4.14) S 2 n{J{x)) = S 2 n-i{x) ■ a:“' 


v(‘) 


i=0 


with the constraint 
(4.15) 


MU 


i—O _7=0 

We also have the corresponding equations for the action of I. 

p-i 

(4.16) S 2 n{I{x)) = S 2 n+l{x) ■ [] 4 


(fc) 


(4.17) S2„+i(/(cc)) = S2 „(t) • n 

i=0 

where Xi is the i-th coordinate of Cx. 

To make relations more uniform, we introduce a slight change of notation: define 
the sequences and {vl^} with the convention that 

(4-18) 0‘2k+l = '^2k+lt ^2k = 

(4-19) P2k+l — '*^2fc+l> /^2fc = ^2k- 

At step n we have 2p -I- 1 variables ((i„, and u)j). 

A first equation simply expresses the factorization: 

p-i 

(4.20) dn= {p- l)dn-i - 

Another set of equations is obtained by expressing that both / and J are involutions: 

(4.21) (p-2)d„ = ^<, i = Q...p-l. 

3¥=i 

It is easy to get from equations (14.2011 and 14.211 

p-i 

< = (7'-2)(in-i-b<_i-^<_i, 

3=0 


(4.22) 


7 = 0 . . .p — 1 . 


8 


J. CH. ANGLES D’AURIACf, J. M. MAILLARDf AND C. M. VIALLET* 


4.5. Singularity structure. We need p additional equations to complete the pre¬ 
vious system. They will be given, under some constraints, by the analysis of the 
singularity structure. The basic idea is that the numbers /3^ (or equivalently 
njj) have a geometrical meaning: they are the multiplicity of some specific 
points on the surface Sn- 

The singularity structure of J is very simple. A singular point is a point whose 
image is undetermined: this happens when all polynomial expressions giving the 
image 03 vanish simultaneously. Any point with more than two vanishing coor¬ 
dinates is singular for J. 

We will look at the action of the pair /, J on the hypersurfaces composing the 
factor Kj of eq. immi. Those are just the n hyperplanes 11^, fc = 0.. .p — 1 of 
equation 

(4.23) Hfe : {xk = 0}. 

All intersections of these hyperplanes are made out of singular points of J. Some 
points are in a sense maximally singular. They are the intersections of all but one 
of the planes 11^, i.e. all but one of their coordinates vanish. There are p such 
points 

(4.24) Pfc = [0,...,0,1,0,...,0], k = 0...p-l 
with 1 in (/c -I- l)-th position. 

The singularity structure of / is the same as the one of J, up to the linear 
change of coordinates C. There are in particular p distinguished singular points 
Qk, k = 1.. .p of I: 

(4.25) Qk=C-^Pk^ fc = 0...p-l. 

To complete the set of equations 62011 , 6211 ), , we need to explore in some 

more details the singularity structure of the maps. What matters is the interplay 
between / and J. 

The map J sends the hyperplane Ilfc onto the point Pk 114.2411 . The 

subsequent images depend on what q is. 

The situation is tractable when g is a prime number, in which case the subsequent 
images of 11^ always go back to the point Pk after a finite number of steps, actually 
one or three steps. There, we meet a singularity, and the equation of Ilk factorizes. 
We will examine the case where q is a prime number, q = 2p — 1. 

The coordinate xq plays a special role and the point Pq behaves differently from 
the other points Pg, s = 1.. .p — 1 . 

Whatever q, the transformation of the hyperplane IIo reads : 

(4.26) IIo ^ Pq —^-fb IIo 

We use the following convention concerning the arrows: when a variety is sent by 
the map onto a variety of same codimension we use the plain arrow —*■. When the 
codimension of the image is lower (blow-down) we use the symbol and when 
it is larger (blow-up) we use the squiggly arrow A blow-up for the birational 
mapping K corresponds to a blow-down for its inverse K~^. 

The action of / and J on the hyperplane IIs reads : 

(4.27) Ugf^Pg^Rg^Rg^Pg^Ug 

The points Rg have coordinates [±I, ±1,..., ±1]. For example for q — 5, R 2 = 
[-I-I, -l-I, —I] and i ?3 = [-I-I, —I, -l-I], while, for q = 7, R 2 = [-fl, —I, —I, -|-I], R 3 = 
[—I, —I, -|-I, -l-I] and R 4 = [—I,-|-I, —I,-|-I]. 
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n 

dn 

IP 


0 

1 

0 

0 

1 

p- 1 

0 

0 

2 

(p- 1)^ 

P-2 

0 

3 

p^ — 3p^ -1- 2p -1-1 

1 

1 

to 

0 

4 

(p - l)(p^ - 3p^ -1-p -1- 3) 

(p- l)^(p- 2) 

P-2 


Table 2. The initial values of dn, and for 0 < n < 4. The 
values for n = 4 have been deduced from the three previous lines 


The pattern is similar for the points Qk ■ It is obtained from the previous one by 
the linear change of coordinates defined by C. The planes Ilfc are replaced by the 
planes 11'^, = C~^Tik and the points Rs are replaced by the points = C~^Rs. 

When q is not a prime number, the pattern is different: the successive actions of 
/ and J leads to singular points other than the P^’s and Qfc’s. In annendixfXIthe 
case g = 9 is studied as an example. 

Relations allow to relate the multiplicities of the singular points 

Pfc on different surfaces Sn- Since Pq ^ Pq in (Eq. 14.2611 we have: 

(4.28) < = 
and since Pg Ps in (Eq. 14.2711 we get: 

(4.29) < = <_ 3 , s = l...p-l 

4.6. End of the proof. The previous analysis shows that when q is prime, the 
factors Xi (resp. W) 1 < i < p appear with the same exponent. In other words, 
for q a prime number, the points Pi, P 2 ... Pp-i play an equivalent role, they will 
have the same multiplicities on each Sn, and we will use Un to denote their common 
value. 

Using II4.28II . 114.2911 together with 114.2011 and II4.22II we get 
dn = (p-l)d„_i-u°_i - (p-l)u),_i, 

(4.30) U°n = {p-2)dn-2- {P-I)ul_ 2 , 

Un = (p- 2) dn-4 - U°_4 “ (P “ 2) U^_4. 

The rate of growth of the d„’s is the inverse of the modulus of the smallest eigen¬ 
values of the 12 X 12 matrix given by the above linear system. The outcome is that 
the complexity of K is the inverse of the smaller root of 

(4.31) x^ + (^2-(p-lf^ x + l = 0. 

To get the full expression of the generating functions, we need to specify the 
initial values of d„, and Un- They can easily be calculated with the help of 
formal calculation software. The results are summarized in table |2l 

Note that when q is not prime, we may still write a set of recursions similar 
to Il4.d0ll . The system is not complete, and cannot be obtained from the analysis 
presented in sectional (see Appendix EJ. 

5. Arithmetical approach: complexity through number oe digits 

The third approach consists in calculating the image of integer points, and eval¬ 
uating the growth of the size of the coordinates, through the number of digits. 
It means that we do not try to calculate the iterates formally. This method was 
already experimented in [^. 
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Obviously the integer coordinates become extremely large, as large as 10™°° 
and we used the library GMP to implement the program |T^. At each step of 
the calculation we factor out the greatest common divisor of the components. We 
assume that the existence of a common factor between all the coordinates is due 
to a factorization of the underlying polynomial. This assumption is valid, at least 
after the first step where accidental factorization could occur. The degree of the 
polynomial is estimated as the number of bits used to store a typical entry (i.e. 
log 2 (My ) ). The algorithm proceeds as follows: i) construct a random matrix of 
integers respecting the equalities of the pattern under consideration, ii) replace each 
term by its cofactor, iii) divide every terms by the greatest common factor of all of 
them, iv) replace each term by the product of all others, v) divide every terms by 
the greatest common factor of all of them, vi) record the number of digits used to 
store the matrix elements. Note that one can exchange ii) and iv) without altering 
the results. The procedure is iterated for as many steps as possible, and possibly 
several runs with different initial matrices are performed. Note that for pattern 
involving only very few variables it can be efficient to write directly the recursion 
relation over the variables. 

The results are summarized in the Table 0 giving the value of the complexity for 
various values of q and for the four patterns introduced above. For cyclic matrices 
and general q it has been shown in ref (3 that the value of the complexity A of 
K = I o J is a. quadratic integer which is the inverse of the smaller root of 

(5.1) + {2-{q-2)^)x + l. 

In Table 0 an empty cell means that we have not been able to compute the 
corresponding A. This is due to the fast growth of the coordinates, preventing us to 
perform a sufficient number of numerical iterations. The number of digits displayed 
is just an indication of the estimated accuracy of our numerical result. When the 
values are known analytically we display six digits. 

6. Conclusion 

The three different approaches we have used give corroborating results. This 
gives us very good confidence both in the heuristic method of section 0 and the 
more numerical approach of section 0 thanks to the proof given in section 0 We 
see by comparing the two last columns of table 0 that Xq happens to be extremely 
close to As, as well as to Ac. This allows us to state the conjecture: 

Conjecture. The complexity of the transformation K = I o J for the general ma¬ 
trices (pattern (G )), for symmetric matrices (pattern (S)), and for cyclic matrices 
(pattern (G) ), are the same. Their common value is the inverse of the modulus of 
the smaller root of — 4 q + 2) a; + 1 = 0. 

Such a result would mean that although the number of parameters of pattern (G) 
and (S) is much bigger than the one of pattern (C), the latter captures the entirety 
of the complexity of the product of inversions K = I o J . This might be related 
to the structure of bialgebra of the set of square matrices equipped with ordinary 
matrix product and Hadamard product. Phrased differently, the “squeleton” formed 
by the cyclic matrices encodes the structure of the whole bialgebra. This deserves 
further investigations which are beyond the scope of this paper. 

Appendix A. The cyclic symmetric case eor g = 9 

We consider in this appendix the case g = 9. Since g is not a prime number, our 
result of section 0 does not apply. 
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q 

Cyclic Symmetric 
Acs 

Cyclic 

Ac 

Symmetric 

As 

General 

Ac 

g Numerical 
Analytical 

1.00026 

1 

6.85424 [7] 
6.854102 

6.85972 [7] 

6.85848 [6] 

g Numerical 
Analytical 

4.0003 [10] 

4 

13.9288 [5] 
13.928203 

13.8811 [5] 

13.965 [4] 

^ Numerical 
Analytical 

6.8541 [7] 
6.854102 

22.9583 [4] 
22.956439 

22.9771 [4] 

22.972 [4] 

g Numerical 
Analytical 

10.3317 [6] 

33.972 [4] 
33.970562 

33.970 [3] 

34.118 [3] 

g Numerical 
Analytical 

12.8326 [5] 

47.027 [3] 
46.978714 

47.040 [3] 

47.000 [3] 

Numerical 

Analytical 

17.9453 [4] 

62.004 [3] 
61.983868 

62.091 [3] 

62.085 [2] 

Numerical 

Analytical 

22.9562 [4] 
22.956439 

79.02 [3] 
78.987340 

79.133 [2] 

80.711 [2] 

^2 Numerical 
Analytical 

25.8105 [4] 

98.03 [3] 
97.989795 

99.17 [2] 

100.32 [2] 

Numerical 

Analytical 

33.972 [3] 
33.970562 

130.3 [3] 
118.9916 

121.6 [2] 

121.5 [1] 

Numerical 

Analytical 

39.169 [2] 

142.8 [2]) 
141.9930 

144.5 [2] 

144.2 [1] 

Numerical 

Analytical 

42.19 [2] 

167. [2] 
166.9940 

170. [2] 


Numerical 

Analytical 

49.10 [2] 

194. [2] 
193.9948 



Numerical 

Analytical 

61.66 [2] 
61.983868 

224. [2] 
222.9955 




Table 3. Complexities oi K = I o J for various values of q, for 
patterns (CS), (C), (S) and (G). The numerical and analytical 
results are displayed. The number in brackets is the number of 
iterations of K we have been able to calculate. 


The number of homogeneous variables is p = {q + l)/2 = ^. We use the same 
notation as in the text for the hyperplane 11^ and the point Pi. In addition we 
define the three points Qi = (1,1, —1, —1,1), Q2 = (1,1,1, —1, —1) and Q4 = 
(1, —1,1, —1,1). We also introduce the codimension-two variety Ilo^a defined by the 
two equations xq = 0 and X 3 = 0. The singularity structure is: 


Ho 

j 

^ To ■ 

I 

>Po 

Ilg 




Hi 


I 

*■ Qi 

^Qi- 

I 

>Pl 

-w III 

TI2 

P2 ' 

I 

>■ Q2 


I 

>P2 

-w II2 

Ha 

J 

I 

*■ Qs 

Ho,3 




n4 

J 

^P 4 - 

I 

> Q4 

^Q 4 - 

I 

>P 4 

4n4 


the subsequent iterates of IIo^s are non singular. We see that there will be six sets 
of exponents, and related to xq, and related to xi, X 2 and X 4 , and 
finally and related to X 3 . The equations expressing the degree drop due to 
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n 

dn 


Wn 

Un 

0 

1 

0 

0 

0 

1 

4 

0 

0 

0 

2 

16 

3 

0 

2 

3 

59 

12 

0 

8 

4 

216 

46 

3 

32 


Table 4. The initial values of dn, and for 0 < n < 4. 


the factorization, and the fact that / and J are involutions, are: 


■^n +1 


^n +1 


^n +1 


^n +1 


= Adn - - ul 

— 3dyj 

= 3d„ - - ul 

= 3dn — ul — Sul. 


Moreover the singularity structure shown above yield : 


■^n +1 


■^n+l 


^n -2 


It is clear that an equation is missing to close the system: 

dn+l = Adn — ul — 3ul — ul 




-^n+l 


— 3dn—l 3u^_i '^n—1 

= 3dn—3 ~ 'U-n—3 ~ ‘^'^n—3 ~ '^n—3 


If we suppose that there exists a recursion relation of the form: 

= a dn—q + b ul_„ + c ul_„ + d ul_„ + e, 


^n+l 


where the a, b, c, d, e, as well as the shift q, are integer constants. The hypothesis 
< 7=1 yields: 


(A.l) 


^n+l 


— 2dyj_i 3 'u^_2. 


Introducing, with obvious notations, the generating functions 
d{s) = 
one easily finds: 


Ui{s) = ^uls^, i = 1, 2, 3 


(4-s^-s®)s (2 - 3) (1 + s^) 

= 1 + ^ Ms) = — 


Ui{s) = 


Pis) 

(2s2_3) s4 


M 2 (s) = 


Pis) 

{3s*-2s‘^- 2 ) s 2 


with 


from which 

fM) = 


Pis) ’ ^ Pis) 

Pis) = (1 - s) • (1 - 3 s - 2 _ s3 2 + 2 / - s®) 

(1 + u + 3 m^ — 3m^)^ 


(1 — u) (1 — 13 M + 2 + u3 + 12 m4 _ 8 y5 _j_ 
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